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^D ' Abstract 

O ' 

fv^ We discuss some special classes of canonical transformations of the extended phase 

space, which relate integrable systems with a common Lagrangian submanifold. Various 

parametric forms of trajectories are associated with different integrals of motion, Lax 

equations, separated variables and action-angles variables. In this review we will discuss 

namely these induced transformations instead of the various parametric form of the 

l/~) ■ geometric objects. 
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Introduction 



Let us begin with a simple example. Consider an ellipse defined by the standard implicit 
J> ' equation 

^" ■ — + yl = 1 

^t ■ a 2 b 2 

O ' 

qv ' One can represent this ellipse by the parametric equations 

o" 

O ■ x = a sin (*) > y = b cos(t) , t € [0, 2-k] . (1.1) 



It is known, there are infinitely many parameterisations of a given curve. For instance, we 
can reparameterise an ellipse by using another parameter 

t= {a 2 + b 2 )t + (a 2 -b 2 ) sin(t). (1.2) 



Construction of different polynomial, rational and other parameterisations of the plane 
jJJ ■ curves is a subject of classical algebraic geometry. 

In classical mechanics the same ellipse may be identified with integral trajectories of 
various integrable systems on the common phase space. In this case the parameter t is the 
time variable conjugated to some Hamilton function H . As an example, the first paramet- 



ric form of the el lipse (1.1) is related to the two-dimensional oscillator, while the second 
parameterisation (|l.2| ) may be associated to the Kepler model. 

In this and several other unexpected situations in mathematics, dynamics is occasionally 
invading mathematical objects in which time is not present in the definition and yet the 
object can be endowed with various dynamical system structures, continuous or discrete. 

Thus, there is a problem of finding suitable parameterisation of curves, which are asso- 
ciated with various integrable systems. What do we know about this problem? 

We consider a mechanical system defined by some Lagrangian function C{q, q) or a 
Hamilton function H(p, q) on the 2n-dimensional phase space M with local coordinates 



{<7jjPj}? = i- According to Maupertuis' principle the extremals qj = 7j(i) of the action 
functional 

5= f C(q(t),ffl)dt, q = ( qi ...,q n ) (1.3) 



I A 

coincide with the extremals of the reduced action functional Sq on the fixed energy surface 



= {H(p,q) = E). (1.4) 

Recall that the Maupertuis' principle was at first enunciated in 1744 M, pj|,p!o|. The modern 
interpretation of the Maupertuis' principle may be found in ^ |2^, j|] . The reduced action 
Sq = J pdq is independent of any evolution parameter. Moreover, even its initial parameter 
t and the corresponding Hamilton function cannot be restored from the reduced action 
problem 0]. Nevertheless, solutions of the corresponding variational problem are the initial 
trajectories qj = "fj(t) in the common nonparametric form 0. For instance, trajectories of 
the Kepler system are conic sections ar^ 1 = 1 + 6 cos<^, which may be ellipses, hyperbolas 
or parabolas at E < 0, E > and E = 0, respectively. 

So, any integral trajectory may be parameterised by using another parameter, t,such 
that 

<fc=7i(*) =&(*), t€[A,B], t£[C,D]. 

According to Maupertuis' principle all the initial trajectories have one nonparametric form 
on the surface Q 2 ™ -1 . Therefore, for all these trajectories we could introduce common 
parametric form qj = f3j(t) and a new local Hamilton function H defined on Q . Two 
main problems are how to find the new parametric form of the initial trajectories and how 
to obtain the new Hamilton function defined on the whole phase space M. 
Consider the corresponding Hamilton-Jacobi equation 

dS H (dS_ dS\ _ Q 

dt \ dqi ' " ' ' dq n J 

In the invariant geometric Hamilton-Jacobi theory |44|] any hyperplane Q 2 " -1 in M. may 
be called the Hamilton-Jacobi equation. Its solution is an 7i-dimensional Lagrangian sub- 
manifold C^") in A4, such that C^ n > C Q 2n ~ 1 . By definition, a Lagrangian submanifold is 
one where the symplectic form $7 vanishes when restricted to it, i.e. fi| c = 0. This def- 
inition is completely invariant with respect to change of local coordinates and parametric 
representations of the Lagrangian submanifold [fl H. As above, we could consider various 
parameterisations of a given Lagrangian submanifold. Each new parametric form yields a 
new Hamiltonian system related to the same geometric surface. 

Below we consider integrable systems on M. with n integrals in involution. Accord- 
ing to the Liouville theorem for any integrable system the corresponding rt-dimensional 
Lagrangian submanifold depends at least of n-arbitrary constants. So integrability is a 
geometric property and it does not depend on the choice of the parameterisation of the 
Lagrangian surface. Starting with a known Lagrangian surface of some integrable system 
we can try to get new integrable models by using various parametric forms of this surface. 
In this case, we can expect that the initial and resulting integrable systems have a lot of 
common properties. The main problem is how to find different parameterisations and the 
corresponding sets of integrals of motion defined on the whole phase space M.. Generally 
there is no rule for how to proceed. Each case is different. 

Usually the Lagrangian submanifold depends on m > n arbitrary constants. The n 



constants a\, . . . , a n are identified with the values of integrals of motion Ij = oy 44 , while 
the remaining m — n constants a\, . . . , a m _„ are free parameters. In this case we have a 
freedom related to the choice of n integrals of motion from the n + m initial constants 
of motion. This freedom permits us to associate a family of integrable systems with one 
Lagrangian surface. In this review we discuss a special class of different parametric forms 
of a given surface, which is associated with mutual permutations of energy E and arbitrary 



parameter a&. The corresponding reparameterisation of C^ n > we will call the generalized 
Kepler change of time. 

The aim of this paper is to bring together some old and new examples of integrable 
systems related with the various parametric forms of the one Lagrangian surface. The 
passage from a given parameterisation to the another one gives rise to the transformations of 
all the properties of integrable systems, such as integrals of motion, Lax equations, separated 
variables and the action-angles variables. In this paper we discuss namely these induced 
transformations instead of the various parametric form of the geometric objects. 

The initial symplectic form Q, is equal to zero on the Lagrangian submanifold C^ n ' . 
We can suppose, that the space M. may be equipped with another form $7, which is equal 
to zero on the same surface C^ n \ In this case for a given Lagrangian surface time t and 
symplectic structure may be changed simultaneously. Moreover, we can try to embed this 
n-dimensional surface into the various phase spaces. For instance, the n-dimensional Kepler 



problem may be identified with the geodesic flow on the sphere S n |17|, 21 . The Kolosoff 
transformation maps the Kowalevski top into the two-dimensional Stackel system |15[ . The 
same top may be related with the geodesic motion on SO (4) pi or with the Neumann system 
on the sphere S 2 |L1| . Here we will not discuss such composed transformations of time and 
phase space. 

By embedding a Lagrangian submanifold C^ n ' into the infinite-dimensional phase space, 
we can identify C^- n ' with an invariant manifold of some hierarchy of nonlinear evolution 
equations. Such finite dimensional manifolds are invariant with respect to the action of 
all flows of the hierarchy and they are naturally expected to be integrable since all the 
flows of hierarchy commute on them. For instance, given Lagrangian submanifolds may 
be realised via stationary or restricted flows. Here we will not discuss relation of finite- 
dimensional integrable systems with soliton equations and restrict ourselves to finding new 
parameterisations of the known trajectories only. 

Below we will consider many well known mechanical systems, such as the Stackel sys- 
tems, Toda lattices, Henon-Heiles and Holt systems, integrable systems with the quartic 
potential and the Goryachev-Chaplygin top. The results we present are, for the most part, 
not new and we do not provide detailed proofs (these can be found in the papers cited). 
What may be new and interesting is an exposition of canonical transformations of the ex- 
tended phase space as different parametric forms of integrable geometric objects and the 
action of these transformations on the properties of integrable systems. 

2 The Maupertuis-Jacobi transformations 

Let M be an n-dimensional Riemannian manifold with the metric gij. On the cotangent 
bundle M. = T*M consider a Hamiltonian system with the natural Hamilton function 
H(p,q) 

n 

H{p,q) = T(p,q) + V(q) = £ 9ij (q) Pl Pj + V(q) . (2.1) 

hj 

On the smooth submanifold Q 2n ~ 1 ( |1,4| ), integral trajectories of the Hamiltonian vector field 
£ = sgrad H (p, q) coincide with integral trajectories of another vector field £ = sgrad H (p, q), 
where the new Hamilton function is given by 

H(p,q) = f(p,q) = J2T^nq~) PiPj ■ (2 - 2) 

i,j 

Integral trajectories have two different parametric forms qj = Jj(t) = Pj(t) on the surface 
Q2n-i on ]y However, the resulting Hamilton function describes geodesic motion and, there- 
fore, on the whole phase space we can determine the so-called Maupertuis transformation 



which relates the initial hamiltonian vector field ^onM with the other hamiltonian vector 
field £ defined on the same phase space M || . 

If t be the time along trajectories of the new vector field £, then the Maupertuis mapping 
gives rise to so-called Jacobi transformations [10[ [H], gTJ of the Hamilton function (2.2) and 
of the time variable 



dt=(E-V{q)\dt. (2.3) 

This Jacobi transformation explicitly describes new parameterisation of the common integral 



trajectories q(t) = q(t) and determines the new Hamilton function H(p,q) (2.2). 

The Maupertuis transformation maps any integrable system with a natural Hamilton 
function H(p, q) into the other integrable system on the same phase space M. Namely this 
property has been used for the search of the new integrable systems (see references within 
Jl9| , no, 0, G71 E8fl). The Maupertuis principle for integrable systems with a nonnatural 
Hamilton function is discussed in [E2| . 

The property of integrability is independent on the choice of parametric form of trajec- 
tories. However, some criteria of integrability drastically depend upon the parameterisation. 
For instance, the method of singularity analysis associates integrability with the Kowalevski- 
Painleve property i.e. the only singularities of the solutions of the equations of motion are 
movable poles (t — to)~ m in the complex i-plane [fy, E3|. There exist cases of integrable 
systems with the rational integrals of motion, whose analytic structure permits solutions 
with algebraic singularities of the type (t — £o)~ m / fc , (k being a positive integer larger than 
one). These systems satisfy to the so-called "weak" Painleve property (see review |23[|). 

We have some examples of integrable systems related with one Lagrangian submanifold 
p3L B8L p9| , which satisfy to the usual Kowalevski-Painleve property and the "weak" Painleve 
property, respectively. These systems are related with the different parametrisations of one 
geometric object. So we can see that a change of parametric form leads to a change of the 
Kowalevski-Painleve criteria of integrability. 

3 Canonical transformation of the extended 
phase space 

To find new parameterisation of the known integral trajectories or of the Lagrangian surfaces 
one has to introduce a new parameter t and the corresponding Hamilton function H . Thus, 
to describe explicitly the following mapping 

t~t, H(p,q)^H(p,q) (3.1) 

we extend initial phase space M by adding to it the new coordinate q n +i = t with the 
corresponding momentum p n +i = —H. The resulting 2n + 2-dimensional space Me [fsl M 
is the so-called extended phase space of the hamiltonian system. We emphasize that H(p, q) 
is the Hamilton function on M, but H is an independent variable in the space Me- The 
energy E is a fixed value of the variable H or the function H(p, q). 

To describe evolution on the extended phase space Me we introduce the generalized 
Hamilton function |lq, p3 

H(pi,...,p n +i;qi ...,q n +i) = H(p,q) - H . (3.2) 

The Hamilton equations for the variables q n +i = t and p n +i = —H are 

dt _ dH 

dr dr 

Here r is a generalized time (parameter) associated to the generalized Hamilton function TL. 
The time variable t is a cyclic coordinate and the conjugated momentum is a constant of 



motion. Other 2n equations coincide with the initial Hamilton equations on the zero- valued 
energy surface 

H(p,q)=H(p,q)-H = 0. (3.3) 

Thus, our initial hamiltonian system on M may be immersed into the hamiltonian system 
on Me- Using this immersion and canonical transformations of the extended phase space 
Me PI p4| we introduce transformations 

(iJ, t) canonical transformations (H , t) 

M E — ► Me 

T I 

M M 

T ~ l 

H(p,q) H[p,q) 

which map an initial Hamilton function H{p 1 q) t—t H{p 1 q) into another Hamilton function 
defined on the same phase space. Of course the similar mapping t h- * t permits us to describe 
new parameterisation of the corresponding integral trajectories. 

Note two different classical definitions of the canonical transformations are known M . 
1. Canonical transformations preserve the canonical form of the Hamilton- Jacobi equations. 
2. Canonical transformations preserve the differential 2-form, Vt — Ej=i dpi A dqt, on M. 
For instance the first definition is used in textbooks on variational principles of classical 
mechanics JlCJ, [16| pj,j6|]- The second definition was later introduced to consider geometry 
of the phase space [gf gf [44| . 

Below we use the first definition of canonical transformations because the Maupertuis- 
Jacobi transformation and the Kepler change of the time preserve the canonical form of 
the Hamilton-Jacobi equation, but it retains the corresponding differential 2-form Q = 
Eri 1 dpi A dq, on the level Q 2 "" 1 (|l|) only @. 

We introduce general canonical transformations of the extended phase space Me 



(3.4) 



1 1— v t , dt = v(p, q) dt , 

H^H, H = v(p,q)- 1 H . 

which change the initial equations of motion 

dqi _i< ,(dq t Tjdv\ dpi _ x ( dpi ~dv 

but preserve the canonical form of the Hamilton-Jacobi equation 

— +H = 0, where S = j (pdq - H dt\ . (3.5) 



and retain the corresponding zero-energy surfaces ( |3.3|) at v(p, q) ^ 

H{p,q) = v{p,q)- 1 U{p,q) = Q. 

Zeroes of the function v(j>, q) determine the behavior of the system with respect to the 
inversion of time pi. Here we do not consider this problem in detail. 



The Maupertuis- Jacobi mapping (2.2), (2.3) may be rewritten as such a canonical trans- 



formation (3.4) of the extended phase space 

T(p, q) h-> T{jp, q) = v(p, qy 1 T(p, q) , v(p, q) = E - V(q) , 

which maps an initial geodesic flow into another geodesic flow. This map preserves integra- 
bility, if the function v(jo,q) is constructed by any potential V(q) 7 which may be added to 
the initial kinetic energy T{p,q) without loss of integrability pi. 



In contrast with the Maupertuis-Jacobi transformations, even if the general canonical 
transformation of M.e fl3.4| ) retains integrability, we have no general method to construct 
new integrals of motion starting with initial ones. To solve this problem in p8|, |39], O, E3, E3 , 
we used some analogies with the one known example of such transformations due to Kepler. 

4 The Kepler change of the time 

We begin with brief description of the Kepler change of the time |l7J . We commence with 
two-dimensional oscillator defined by the Hamilton function 

H osc (p,q)=pj+p1 + a(<£ + qj) + b, a,beR. 



For this system the Kepler canonical transformation (3.4) of M.e with the function 

v(p, q) = q\ + q\ (4.1) 



preserves integrability. After change of the time (3.4) and the point canonical transformation 
to other variables 

x = qiq2 , V = {q\ - q 2 2 )/2 (4.2) 

integral trajectories of the oscillator come to be trajectories of the Kepler problem defined 

by 

H kepl (p,x) = g^j) =p £+pg+ b +,. (4.3) 

Various parametric and nonparametric forms of the common trajectories are discussed in fL7| , 
@, B|. Coincidence of the integral trajectories may be regarded as a local result, whereas the 
corresponding canonical transformation of the extended phase space preserves integrability 
in the whole initial phase space. 



As for the Maupertuis-Jacobi transformation, the function v(p,q) (4.1) in the Kepler 



transformation (4.3) could be identified with the oscillator potential V{q). Below we prove 



that it is a simple coincidence. Nevertheless canonical transformations of the type (3.4) 
have been called the coupling constant metamorphoses in |12, E3] because of this casual 
coincidence. 

The Kepler change of time has been generalized by Liouville |l8f . The Liouville inte- 
grable systems are systems with the natural Hamilton function 

H(p,q)=T+V, 

where the kinetic and potential energies are given by 

n n n 



T = v- 1 {q)Y d <MPh V = v- x {q)Y,Ui, v(q) 






=i 



The functions a%, Vi and Ui depends only on the variable q h . 
For the Liouville system the following quantities 

I j = OjPj + Uj-Hv j , j = 1, . . . , n , 

are integrals of motion in involution and ^ Ij = 0. Thus, we have n quadratic integrals 
of motion including the Hamilton function and consequently the Liouville systems are com- 
pletely integrable. Note the two-dimensional oscillator and the Kepler model belong to the 
Liouville family of integrable systems. 

We recall how equations of motion were integrated in quadratures by Liouville |lq ]. 
From Ij — ctj , one obtains a system of differential equations 

dqj dt 

j = l,...,n. 



sJajicij + Evj-Uj) 



v(q) ' 



Here H = E and time t is associated with the Hamilton function H. Choosing a new time 
variable t according to 

dt = v-\q)t 

we come down to the system of equations 



\/ a A a J+ Ev 3 ~ U j) 



dt. 



It allows one to find integral trajectories qj = "fj{t). After that the new parameter t may be 
expressed in terms of the initial time t by the quadrature 

t = J w(«i(r),...,? n (r))dT. (4.4) 

This transformation is related to a new parametric form of the same trajectories qj = 0j (t) . 
In fact, Liouville tacitly used canonical transformation of the extended phase space 



(3.4) and considered a new integrable Hamilton function H — v(q) H instead of the initial 
one. After integration of equations of motion for the new system we have change parametric 
form of the trajectories in order to describe solution of the initial problem. For instance in 



parametric form of the ellipse the Kepler time (1.2) is explicitly the Liouville quadrature 



(4.4) 



Note canonical transformations of the extended phase space (3.4) have a natural coun- 
terpart in quantum mechanics. Namely, it is known that the standard eigenvalue problem 
of the Hamiltonian operator H (p, q), 

H(p, q)V = {H Q + aV + b)y = E^ , 

may be associated with the eigenvalue problem of the charge operator a 

H(p,q)V = {H Q + {b - E)V~ 1 )^ = -a*. 

In quantum mechanics such a duality of the two eigenvalue problems has been used by 
Schrodinger and many other |25| . Canonical transformation of the extended phase space 



(3.4) is an analogue of this duality. It is interesting that for the first time this duality 
has been studied for the quantum Kepler problem. In the Birman-Schwinger formalism 
function v(q) is called a "sandwich" potential [E5|. Recall that in the Birman-Schwinger 
formalism we can estimate the spectrum and eigenfunctions of the one Hamiltonian, H, by 
using the known spectrum and eigenfunctions of the dual Hamiltonian, H . Moreover, for 
some quantum models it is a single known way to find solutions of the initial Schrodinger 
equation. Below we briefly discuss a similar property for the quantum Toda lattice. 

5 The generalized Kepler change of time 

The Maupertuis-Jacobi mapping is traditionally used for the search of new integrable sys- 
tems. The Kepler change of time and the Liouville reparameterisation of trajectories have 
been used for integration of equations of motion. In this Section we propose some generali- 
sations of the Kepler-Liouville results, which may be useful for the search of new integrable 
systems as well. 

All the Liouville systems are particular case of Stackel integrable systems. Therefore let 
us briefly recall some necessary facts about Stackel systems |33| . The nondegenerate n x n 
Stackel matrix S, its j'-th column of which, Skj, depends only on qj 

dctS^O, ^=0, j/m 
oq m 



defines the set of functionally independent integrals of motion, {/&}£_!, where 

J fe = J2 c i k ( p2 J + U i &') ) ' c J' fe = d^s ' ^'^ 

i=i 

which are quadratic in the momenta. Here C = [cjk] denotes the inverse matrix of S and 
S fcj ' is the cofactor of the element Skj ■ 



Proposition 1 M If the two Stdckel matrices S and S be distinguished by the m-th row 
only, i.e. 

Skj = s^ , k ^ m , 



i/ie corresponding Hamilton functions I m and I m (5.1) with a common set of potentials U. 



j 



are related by canonical transformations of the extended phase space M.e 

im >-> I m = v^ 1 (q) I m (p, q) , dt m = v(q) dt m , (5.2) 

where _ 

detS(gi,...,g n ) , , 

u(qi,...,g„) = — -— -. (5.3) 

detS(gi,. ..,q n ) 

The proposed generalization of the Kepler transformation maps one integrable Stackel sys- 
tem into another integrable Stackel system. For instance, the Stackel matrices for the 
oscillator and the Kepler problem are 

\\), S=(f *). (5.4) 



It is obvious that the Kepler change of time (4.3) coincides with the proposed mapping (|5.2|) 



We return to the Kepler transformation (4.3) of M.e- After permutation of coordinates 



and momenta (51,2 * r ~* Pi,2) the Hamilton function for the Kepler problem 
~ _ p\ +pl + a- l {ql+ql) + a- l b _ H osc 

-^kepl — & 9,9 ^ tt 

PI + P 2 H free 

becomes a ratio of the Hamilton function, H osc , for the oscillator and the Hamilton function, 
Hf ree , for the free motion. 

So for any two integrable systems the ratio of their Hamilton functions could be the 
Hamilton function of a third integrable system on the same phase space. The main remaining 
problem is a search of a complete set of integrals of motion. 

We consider two integrable hamiltonian systems on the common phase space M. . These 
systems are defined by the two sets of independent integrals of motion, {/,-}" =1 , and { Jj}" =1 , 
in involution, i.e. 

{Ij,I k } = and {Jj, J fc } = 0, j, k = 1, . . .,n. 

Introduce the antisymmetric matrix K. = (I <S> J), which is the inner product of the two 
independent vectors of integrals I and J in W l . Any column or row of this matrix defines a 
set of n — 1 independent functions 

Kij = (I (g) J)y = Ii Jj - I } ■ Ji , i, j = 1, . . . , n . 

Proposition 2 H / If all the differences of integrals of motion (Ij — Jj) with the common 
index j = 1, . . . ,n are in involution, i.e. 

[l j -J j ,I k -J k )=Q, j,k = l,...,n, (5.5) 



then the ratio of integrals 

K m = ^ (5.6) 

and n — 1 functions Kj, j ^ m 



"m "T? 



K j = ^-= ±m ° 3 1 l3 " m =K m J i -I i , m*j = l,...,n (5.7) 



are integrals of motion for new integrable system on the same phase space. 



Thus the mapping (5.6) defines a canonical transformation (3.4) of the extended phase space, 



which preserves the property of integrability. To apply this transformation we have to find 



two integrable systems satisfying condition (|5.5| ) 



We consider a pair of the Stackel systems with a common Stackel matrix S and with 



different potentials Uj. Namely, in addition to the system with integrals {Ik} (5.1), we 
introduce the second integrable system with the similar integrals of motion, 

n 

Jk = J2 c i k ( p2 J + W i fe) ) . k = l,...,n. (5.8) 

At least one potential Uj (qj) has to be functionally independent of the corresponding po- 
tential Wj (qj). 

Proposition 3 / pj| / Any two integrable systems defined by the same Stackel matrix S and 
by f unctionally independent potentials Uj (qj) and Wj (qj) satisfy the necessary condition 



(5.1) of the previous proposition. Thus, the ratio of the two Stackel integrable Hamiltonians 



defines new integrable system 

(I m ,Jm) !-* K m =v(p, q)' 1 I m = -^- ■ (5.9) 

It is obvious that all the integrals Ik and Jk differe by the potential part 

2 

(h - Ju) = Y. c ^WAi 3 ) -w 3 (ii)] 

depending on the coordinates q only. Thus systems with a common Stackel matrix S satisfy 



condition ( |5.5| ) 



The Hamilton function (|5.9| ) has the following form 

fj( \ V E"=lCjm[p] + Uj(q 3 )] -An 

H[p, q)=K m = -^—i , /3 m G R . (5.10) 

Ej=l c jm [pj + Wj (qj)\ 

This Hamiltonian H (p, q) is a rational function in the momenta, but next one can try to 
use canonical transformations to simplify it. Occasionally, one obtains again a natural type 
of Hamilton function. For instance, according to pif , integrable systems with the following 
Hamilton functions 



Hi = p k x p k y +a{xy) k + l , a,ke 



k_ 

En = p k x + P k y +a(xy) k + l 



(5.11) 



belong to the proposed family of generalized Stackel systems. At fc = 1 the first Hamiltonian 
coincides with the Hamiltonian of the Kepler problem. At k = 2 the second integrable 
Hamiltonian has been found by Fokas and Lagerstrom (see references within |2J|). In this 
case both initial systems satisfy to Kowalevski-Painleve criterion, whereas the resulting 
Fokas-Lagerstrom system admits asymptotic solutions with fractional powers in t [p3[ . 



6 Properties of the change of time for the Stackel sys- 
tems 



For the Stackel family of integrable syste ms we proposed two different examples, (5.2) and 



(5.9), of the canonical transformations (3.4) of the extended phase space Me- Now we 
consider some properties of these transformations. 

Recall that the common level surface of the Stackel integrals Ij 

M a = {ze R 2 " : U{z) =on, i = l,...,n} (6.1) 

is diffeomorphic to the n-dimensional real torus. We can immediately construct the one- 
dimensional separated equations 

/op \ 2 n 

P 2 j = f -q^ J = Pj(lj) = ^2 aiS v(lj) ~ u j(lj' a ) ' a k eR. (6.2) 

Here Sq is a reduced action functional |33| . Integral trajectories qj(t, cti, . . . , a n ) are deter- 
mined from the following equations 

§/« = - *-■ 

In fact the polynomial P(A) and the contour of integration depend upon the values, ay, of 
the integrals of motion, which are dropped in the notation. 

For the rational entries of S and rational potentials Uj (qj ) the Riemann surfaces (|67 
are isomorphic to the hyperelliptic curves 

2g+l 

C r . ( M 2 j =P(X j )=J2^^- (6-4) 

fe=i 

Considered together these curves determine the n-dimensional Lagrangian submanifold in 
the phase space M = R 2n 

CW : C 1 { Pl ,q 1 )xC 2 {p 2 ,q 2 )x---xC n ( Pn ,q n ), (6.5) 

which is decomposed into plane curves. Applying Arnold's method |2| we find that the 
action variables have the form 



y/P^i) d\j ■ (6.6) 



The Abel transformation linearizes the equations of motion on the Lagrangian submanifold 
C^ n ' in terms of abelian differentials of the first kind on the corresponding spectral curves 



We consider a pair of the Stackel systems related by the first generalization of the Kepler 



mapping (5.2) such that the potentials Uj(X) — J2 a k^ k depend on arbitrary parameters 



a.fc. According to []38| initial and resulting integrable systems are associated with algebraic 



hyperelliptic curves (3.4) C and C described by 

C: i? = J2 o.j X J + a m X m + E\ k + ^ a t X l , 

(6.7) 
C: i? = Y^ aj A j + EX m + a k X k + ^ a t X z . 

Here the n coefficients {aj} and {5j} are values of the integrals of motion such that E = a\ 
and E = a\. The other coefficients aj s K are arbitrary parameters (charges), which 

10 



define the potential part of the Hamilton function H(p,q) = T(p) + V(q,a). Canonical 
transformations of Me give rise to mutual permutation of the energy E and one of the 
parameters a& (charge). 

The initial and resulting Riemann surfaces are topologically equivalent. One can prove 
that the corresponding integrable systems are topologically equivalent too. Moreover, initial 
curves coincide with the resulting curves at the special values of integrals of motion 



E = a. 



E = a k 



'.i > 



J = 2, 



,,72. 



In this case we have two different parametric forms of the common Lagrangian submanifold 
C^ n ' (6.5) depending on n values of integrals oij and constants a k - On the corresponding 



submanifolds M a and M a (6J) integral trajectories of the initial system, qj(t) (3.3), co- 
incide with integral trajectories of the resulting system, qj(t). In the neighbourhood of the 
intersection of these submanifolds we can introduce the common set of the action-angle vari- 
ables (6.6). In this region the Kepler transformation (5.2) retains the differential 2-form f2 
in Me and the function v(p, q) = v(s) depends on the common action variables only. 

We consider a pair of the Stackel systems related by the second generalization of the Ke- 
pler mapping (5.9). According to Bl| initial and resulting integrable systems are associated 
with algebraic hyperelliptic curves (|6.4|) C and C of the form 



C: 



M 



J2a t X l + J2^X 3 



C 



£fcA*+5> 



Here the resulting coefficients are rational functions of the initial ones [|41j. So canonical 
transformations of Me (5.9) give rise to transformations of the modulus of the curves only. 
As mentioned above the corresponding Riemann surfaces and integrable systems are 
topologically equivalent. At the special choice of integrals {a, 5} and parameters of po- 
tentials {cii, b} integral trajectories of the initial system coincide with trajectories of the 
resulting system. In a small region of M this generalization of the Kepler transformation 



(5.9) preserves the differential 2-form Q in Me and the function v(p,q) depends on the 



action variables only. 

In an attempt to u nder stand the origin of the conservation of integrability by canonical 
transformation of Me (^4) one can attempt to rewrite equations of motion in the Lax form 



{H(p,q),L(\)}= [L(X),A(X) 

For the some classes of Stackel systems the Lax matrices have been constructed in |37], |3^, E(J . 
We consider some examples only. One of the simplest Lax matrix L(X) for the one- 
dimensional Stackel systems is given by 



L(X) 



( 



\ 



V 

4>(X) 



X-q 



MN 



X-q \ 
/ 



-P 



(6.8) 



Here </>(A) = ^ 4>k X k is a parametric function of the spectral parameter A and the elements, 
[z]mn, arc the linear combinations of the Taylor projections 



[z]n = 



£ z k X k 



_fc= — oo 



J N 



N 

E 

k=0 



z k X k , 



(6.9) 



or the Laurent projections |3(|. The coefficients <pk of the function </>(A) arc the constants 
of motion, which may be parameters a k or integrals of motion. So at arbitrary M, N the 
family of the Lax matrices L c j > (X) may be associated with a web of algebraic curves instead 
of one concrete curve. 
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For instance, we present some one-dimensional Hamilton functions, the functions c^(A) 
and the corresponding Lax matrices associated with the first Taylor projection [z]qi- An 
application of the pure numeric function 4>(X) yields 



H = p 2 



aq 



bq, 



-aX 2 -bX. 



L = 



P A-g 

—a (A + q) - b -p 



A = 



f 
-a 



(6.10) 



The corresponding spectral curve is given by 



C 



(i 2 = aX 2 + bX- H . 



Note that in the spectral curve we can substitute integrals of motion, while in the corre- 
sponding Lagrangian submanifold we have to substitute the values of integrals of motion. 
From (6.7) the first possible change of this curve looks like 



C i-> C : /i 2 = aX 2 + {b- H)X + c. 
The associated canonical transformation of the extended phase space, 

H = v- 1 {H + c) =q- 1 (H + c), 
changes the Lax matrices by the following rule 

<t>= -a A 2 - (b-H)X, L = L 

The second possible change of the curve 

C i-> C : p? = (a-H)X 2 + bX + c 
may be related to the other canonical transformation of the extended phase space, 



H 





1 



A = v- 1 (q)A. 



(6.11) 



H = v- 1 (H + c) 



(H + c), 



such that 



{a-H)X 2 -bX. 



(6.12) 



L = L + H 




X + q 



A = v-\q) 



A + H 
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It is not hard to check that the Poisson bracket relations for all these Lax matrices L(X), 
L(X) and L(X) are closed into the linear r-matrix algebra. At the first case the r-matrix 
is a constant matrix r = 11/ (A — /i), whereas the second and third r-matrices depend on 
dynamical variables. Here n is a permutation of auxiliary spaces ||. 

We turn now to the original change of time in Kepler problem (fL3) . The Lax matrices 
for the two-dimensional oscillator, 



£(A) = 



Li(X,pi,qi) 

L 2 (A,P2,<?2) 



A(X) = 




(6.13) 



may be constructed from two independent 2x2 blocks Lj(X) and Aj(X) of ( 6.10 ). The 
corresponding spectral curve C — C\ x C2 is a product of two hyperelliptic curves. 



12 



The Kepler mapping (4.3) gives rise to the following transformation of the Lax matrices 
£(A) = C(\) — Hkepi 



( \ 

A + 9 i 



\ X + q 2 J 



(6.14) 



A(X) = v(p, q y 



-4(A) - H, 



kepi 



^ o o o o \_ 

10 

oooo 
\ o o i o J 



The spectral curve C remains a product of two new hyperelliptic curves. The initial oscillator 
and the resulting Kepler model are separable in the common variables, which lie on these 
curves. Note that these separated variables are cartesian coordinates for the oscillator and 

parabolic coordinates for the Kepler problem. 

Similar transformations of the Lax matrices ( |6.1l| ), ( 6.12 ), ( 6.14 ) may be proposed for 
the other two-dimensional Stackel systems separable in cartesian or parabolic coordinates. 
For the two-dimensional Stackel systems separable in elliptic or polar coordinates transfor- 
mation of the Lax matrices has a more complicated form |38| . These transformation may 
be constructed by using two different outer automorphisms of infinite-dimensional represen- 
tations of underlying sl{2) algebra proposed in Go). 



In the next Sections we consider similar transformations of the Lax matrices (6.11) and 
the spectral curves (6.7) for integrable systems associated with non- hyperelliptic algebraic 
curves and for the non-Stackel integrable systems. 



7 On integrable systems with quartic potential 

According to p3, [39| |43| , canonical transformations of the extended phase space relate three 
integrable cases of Hcnon-Heiles systems with the three integrable cases of the Holt systems. 
One of these systems admits a 2 x 2 Lax matrix and the corresponding spectral curve is 
an hyperelliptic algebraic curve. Another two systems possess 3x3 Lax matrices and the 
corresponding spectral curves are the trigonal algebraic curves /i 3 = P(X). As for the Stackel 
systems transformations of th ese 3 x 3 Lax matrices are shifts of their entries by the element 
of the extended phase space fl6.1l|) . 

We consider two integrable systems with quartic potential for which 4x4 Lax matrices 
were constructed in J4| by applying relations with stationary flows of some known inte- 
grable PDEs. The first system belongs to the Stackel family and is separable in cartesian 
coordinates. Its Hamilton function is 



H(p,q)=pl 



2 1 
p 2 2 + - 4 



'A 



r 2 2 

6 q l q 2 



4\ 



(7.1) 



As above we could construct some block 4x4 Lax matrices (3.13) for this system in cartesian 



separated variables. The corresponding spectral curve is a product of hyperelliptic curves. 
The two pairs of the separated variables lie on these two curves. 

Another 4x4 Lax matrix for the same system may be obtained from the Lax represen- 



ts 



tation for the Hirota-Satsuma coupled KdV system B. This Lax matrix 



/ 



L(X) 



V 



-A 



-piqi 


fl? 





1 \ 


\-p\ 


Pi 9i 


ql+ql 

2 








A 


-P2?2 


ql 


ql + ql 

2 





A -Pi 


P2<?2 y 



does not have a pure block-diagonal structure and the corresponding spectral curve /i 4 = 
P(A) is not a product of two hyperelliptic curves. The second Lax matrix A(X) may be 
found in 0. 

For this Stackel system canonical transformation of the extended phase space ( |5.2|) 

^=(9?+g2 2 ) _1 (^-6) 

gives rise to the shift of the first Lax matrix and the rescaling of the second Lax matrix 



L = L + H 



^ ^ 

10 



\ A J 



A = ir 1 A. 



In contrast with the Stackel systems (6.7), transformation of the corresponding spectral 
curves acts on the both side of the equation defining the curve 



C 
C 



I' 



/" 



= A 3 -HX 2 +JA, 

; -2fi 2 XH =A 3 -{H 2 -b)X 2 +JX. 



(7.2) 



Note that after canonical transformations of the other variables (p, q) the new Hamilton 
function H may be rewritten in the natural form 



H = p 



2 , 2 x z + 2y 

Py ■>■ 



2 y/x 2 + y 2 



The second integrable system with the quartic potential is the non-Stackel system de- 
fined by the following Hamilton function 



H = pI+pI-\ (qi + 6qlq 2 2 



»<£ 



(7.3) 



This system is separable after a so-called quasi-point canonical transformation |2j 
The second system possesses 4x4 Lax matrix 



/ 



L = 



qiq\ , 

2 + <7iPi 


-ql 


+ pf +qiq 2 Pi + 2 


qiq\ 

2 


-2q 2 X 


2 A 


*(-P2 + q 2 2 + q 4) 






q\P\ 



2q 2 





9 ql 
P2 + ql + -j 




qiql 

2 


+ qiPi 


2 2 
A 


+ p\- 


- qiq2Pi 



-«1 



A g 2 <7? 



<7iPi 
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which was obtained using a gauge transformation of the Hirota-Satsuma coupled KdV system 



For this system we propose the canonical transformation of the extended phase space 
given by 

H = (ql+4q 2 2 )- 1 (H-b), 

which preserves integrability. As for the Stackel systems, transformation of the Lax matrices 
retains a shift of the first Lax matrix and rescaling of the second Lax matrix 



L = L + 2H 



/ \ 

10 



\ A ) 



A = %r 1 A. 



The spectral curves are changed by the rule 



C: 

C: 



/i 4 =A 3 +AHX 2 +JX, 

/i 4 -8n 2 \H =A 3 +4(6 -AH 2 ) A 2 +JX. 



(7.4) 



Recall that for Stackel systems with quadratic integrals of motion canonical transformations 
of the extended phase space (5. 2), (5. 9) preserve separated variables lying on the common 
hyperclliptic curve (3.7). 

For the Henon-Heiles systems and systems with quartic potentials transformations of 
the trigonal spectral curve /it 3 = P(X) |3j| and the curve /i 4 = P(X) (7.2, 7.4) have a more 



complicated form. In this case separated variables for the initial and resulting systems are 



different 1 39 . It means that the proposed canonical transformation of the extended phase 
space preserve integrability, but seriously changes other properties of the systems. 



8 The Toda lattices. 

Before proceeding further, it is useful to recall some known facts about the Toda lattices 
(all details may be founded in the review pq). 

Let q be a real, split, simple Lie algebra of rankg = n and let K(-,-) be its Killing 
form and P be a system of simple roots, we identify the phase space with the coadjoint 
algebra M ~ Q* R — Q* + (B g*_. Here g+ is a Borel subalgebra and g_ is the opposite nilpotent 
subalgebra of g. 

The Lax matrices for the Toda lattices are 



^2 PiK + ^2 ' ex P K ( a -. l) 



a£P 



a£P 



A 



^2 expK(a,q) ■ e. 



aEP 



The Hamilton function is given by 



H(p,q) = ±K(£,£) = ±K(p,p) 



a a e 



22 a ae a 



(v) 



aeP 



Recall that in a shifted version of the Adler-Kostant-Symes scheme in order to construct 
the Toda orbit (|8J ) in M we have to translate a dynamical orbit living in q* + by adding to 



it a constant vector e = ^2a a e a from the remaining part of M.. This vector e has to be a 
character and has to be a constant. 

We replace now the phase space M by the extended phase space Me- Roughly speak- 
ing, to consider the Stackel systems we exchange the pure numerical function <fi (|6 .10) by 
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the function with coefficients from M.e (6.11), (6.12). By a similar reasoning we try to 



construct the same modification of the Adler-Kostant-Symes scheme. Namely we translate 
the Toda orbit in M. ~ q* r by adding to it a constant vector from the remaining part of the 
whole space M.e- As above this vector has to be a character and has to be a constant with 
respect to the new time. 

In addition we impose a constraint on the possible change of parametric form of trajec- 



tories. As for the Stackel systems (6.7) the initial invariant polynomial has to generate the 
arbitrary constant 

K{Z,Z) = -b (8.3) 



instead of the Hamilton function (8.2). This condition together with the form of transfor- 



mations of the Lax matrices dictates to us a very special choice of the functions v(p,q) in 



(3.4) for the Toda lattices. 



Proposition 4 M/ For each simple root (3 € P and for any constant fc^et the following 
canonical transformation of the extended phase space M.e 

dt = e P{q) ■ dt , 



H a = e 



-(9(9) 



(8.4) 

maps the Toda lattice into the other integrable system. This canonical transformation induces 
the following transformation of the Lax matrices 

£ p = £-H p -ep, A = e-K*>-A. (8.5) 



Here H , C and A are the Hamiltonian (8.1) and the Lax matrices (8.1) for the corresponding 
Toda lattice. 

In this proposition we explicitly determine the new parameter t and the new associated 
Hamilton function H . The corresponding spectral curves depend on the choice of a repre- 
sentation of q. We prove that the Toda lattice and the new integrable system relate to a 
common geometric object in the one example only. 

The number of the functional independent Hamilton functions Hp, (3 6 P depends 
upon the symmetries of the associated root system. For closed Toda lattices associated with 
the affine Lie algebras canonical time transformation has a similar form. Similar canonical 
transformations may be applied to the relativistic and discrete time Toda lattices. 

We describe explicitly some new integrable systems related to the standard three- 
particle Toda lattice and two-particle Toda lattices associated to the affi ne algebras X2 . 



After an appropriate point transformation of coordinates ( similar to ( |4.2| ), see |4lJ]) all the 
Hamilton functions have the common form 

H = p x p v + — + ax Zl y Z2 +cx Sl y S2 +d, a,b,c,dER. (8.6) 

xy 

where Z\ 2 and S\ 2 are the roots of the different quadratic equations and are related to the 
angles of the corresponding Dynkin diagrams. Below we show these equations explicitly 



B 



A 3 ■ 


z 2 + 3z + 3 = 


s 2 + 3s + 3 = 


(i) r w . 

2 U 2 ■ 


z 2 + 4z + 5 = 
z 2 + 2z + 2 = 


s 2 + 4s + 5 = 
s 2 + 3s + 5/2 = 


D±»: 


z 2 + 2z + 2 = 
z 2 + 2z + 2 = 


s 2 + 2s + 2 = 
(s + 2) 2 = 


G«: 


z 2 + 2z + 4 = s 2 + 5s + 7 = 

z 2 + 2z + 4 = s 2 + 3s + 3 = 

z 2 + 3z + 7/3 = s 2 + 3s + 3 = 
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Originally the integrable system with the Hamilton function H ( |8.6| ) associated to the root 
system A 3 was found by Drach M . 

The corresponding second integrals of motion K are polynomials of the third, fourth and 
sixth order in momenta. Note that for the algebra A 3 all the three Hamiltonians Hp, (3 G P 
are equivalent. Two different Hamilton functions are associated with the algebras B 2 , C 2 



J 2 1 w 2 

iW T^v +lna fit 1 *) 



and D 2 . For the G 2 algebra we have three independent potentials in (B.6) 



9 On the common properties of the Toda lattices and 
the dual systems. 

For the periodic Toda lattice associated with the root system A n the Hamilton function is 

1 ™ 

Here {pi, qi} are canonical variables and the periodicity conventions qi+ n = q% and pi+ n = Pi 
are always assumed for the indices of qi and pi . 

The exact solution of the equations of motion is due to existence of a Lax equation with 
the following nxn Lax matrices |L|, ^ 

n n—1 

i=l i=\ 

+ ne^-^E ltn + a n iT l E nil , (9.2) 

i=\ 

where Ei^ stands for the n x n matrix with unity on the intersection of the ith row and the 
kth column as the only nonzero entry. 



According to J4l|, |42j canonical transformation of the extended phase space (3.4) by 

v(p,q)=exp(q j -q j+1 ), H = e^ 1 ^ (H - b) , b e R (9.3) 

maps the Toda lattice into the dual integrable system with the following equations of motion 

dq t _ lf sdqi dpt _ lf .dpi ~ 

Tf =v (q)-, -^^v (q) — + H(6 i , j -6 i , j+ r). (9.4) 



Associated with the different indexes j canonical mappings ( |9.3| ) are related with each other 
by canonical transformations of the other variables (p, q) . 



The mapping (9.3) gives rise to the following transformation of the Lax matrices 

C(n) = £(n) - ft E jti+1 , A( f i)=v- 1 (q)A(u). (9.5) 

The corresponding transformation of the spectral curves is 



n— 3 



(9.6) 
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Here p = J\ = J2 Pi ^ s the total momentum, H and H are the corresponding Hamilton func- 
tions and Ji, Ji are integrals of motion. Substituting the fixed values of integrals of motion 
in the product of curves one can construct the corresponding Lagrangian submanifold. 

Applying Arnold's method fl pf to the standard form of the hyperelliptic curves C and 
C (|9.6|) one construct the action variables 



1 



P(A) 



P(A) 



\ 



P(A) 2 -4j]a, dX, 



(9.7) 



\ 



P(A) 2 -4(a, -H)Y[ ai ] dX, 



where Ai and Ai are A-cycles of the Jacobi variety of the algebraic curves ([Tq), respectively 
||. In fact the polynomials P(A), P(X) and ^4-cycles depend on the values of constants of 
motion, which are dropped in the notation. The Abel transformation linearizes equations 
of motion in terms of first kind abelian differentials on the corresponding spectral curves. 

Let parameters a^ determine potential of the Toda lattice fl9.l| ) and parameters a,{ and b 
define the potential of the dual system ( |9.3| ) . At the special choice of the values of integrals 
of motion 

IT a* 



H = b. 



H — a, ■ — 



n 



Ji Ji — Oij, 



(9.8) 



i±j 



the initial curve C is equal to the resulting curve C (9.6). Thus, as for the Maupertuis 
-Jacobi mapping B, integral trajectories of the Toda lattice coincide with the trajectories 
of thtMiual system on the intersection of the corresponding common levels of integrals M. a 
and M. a . In the neighbourhood of this intersection we can introduce the common set of the 
action variables ( |9.7| ) for the both systems. In this small subvariety of the phase space the 
function v(p, q) = v(s) depends on the action variables only. 

At etj = 1 the Poisson bracket relations for the n x n Lax matrices can be expressed in 
the r-matrix form 

{£(//) , C(v) } = [ r 12 {n, v) , L{ix)\ + [r 21 (/x, v) , l{v) ] . 



Here we used the standard notations 

C(ji) = C{ii) <8> I , C{v)=I®C{v), 



r 2 \ (Mj v) = -n ri2 (v, n) TL ■ 



and n is the permutation operator in C" x C" [pi . Canonical transformation of the extended 
phase space (9.3) maps the constant r-matrix for the Toda lattice 



r 12 (^)=r-» st (A^) 



jj, — v 



v Yl + v Yl E " 

m>i m<i J 

into the following dynamical r-matrix 

foG*, v) = r^ nst (M, v) + A{v, q) ® E jd+1 , 

where the second Lax matrix A{v, q) and, therefore, the dynamical r-matrix rij(/J,, v) depend 
on coordinates only. 

Another known 2x2 Lax representation M pll for the same Toda lattice is equal to 



T(A)=ii(A)---i„(A), 



dT 
~~dt 



T{X),A n (X) 



(9.9) 
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where 



u = 



X + Pi 
-a,i-i e~ 







A = 



A e«* 

-a.e" 9 - 1 



(9.10) 



Canonical transformation of the extended phase space (^3) gives rise to the following trans- 
formation of the Lax matrices 



T(X) 



L\ ■ ■ ■ Lj_ 



i-i 



Lj Lj+i 



H -b 




Lj+2 ■ ■ • L n , 



A n (X) = v- l {q)A n {\) 



(9.11) 



At at — 1 the Poisson brackets relations for the 2x2 Lax matrices T(X) ( |9.9|) satisfy the 
following Sklyanin r-matrix relation 



{T(X) , T{v) } = [R{X-v), T(u)T(v) ] , R(X - u) 



n 



X-v 



(9.12) 



The mapping (9.3) transforms these quadratic relations into the following polylinear ones 



{r(A),f(i/)> 



1 2 

= \R(X-v),T{X)T(y)\ 

+ [rit n {\v),f{X)] + [4l n {X,u),T^)'. 



The corresponding dynamical r-matrix is given by 

rf% n (X,v) =A n (X,q) ® f L x • ■ -Zy_i ■ ( Q Q J • L j+1 ® L n 

Here all matrices Lk depend on the spectral parameter v and A n {X, q) is the second Lax 
matrix (|9.9|). 

Now we look at the separation of variables in the framework of the traditional consid- 
eration of the Toda lattice. A complete list of references can be found in |13, pi BiJ . Below 
we put Oi=l and j — 1 without loss of generality, so that 



H = expfe - qi ) (H + b) , T 



L X L 2 + 



H + b 





This transformation changes the first row of the Lax matrix T(A) only 



• L-x- ■ ■ L r . 



T = 



A(A) B(A) 
C(A) ©(A) 



/ A(A) 1(A) 
C(A) D(A) 



(9.13) 



The separated variables {Ai A2 . . . , A n _i} for both systems are zeroes of the nondiagonal 
common entry 

n-\ 

C(X)= 7 -l[(X-X i ). (9.14) 

An additional set of variables is defined by the second common entry 

(Ui = D (Aj) such that {Xi,\og[ik} — Sik i,k = 1, . . . ,n— 1 . 

From detT(A) = 1 and detT(A) = (1 — H) one immediately obtains the one-dimensional 
equations 



A(Ai) = ^ 

A(X t ) = (1 - H)^ 1 






(9.15) 
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At the special choice of parameters and values of integrals (9.8) initial separated equations 
coincide with the resulting ones. 

Thus we prove that the initial and resulting systems have a common set of separated 



variables. On the other hand canonical transformation of the extended phase space (9.3) 
changes the form of the Backlund transformation |42| and the bihamiltonian structure, which 
are known for the Toda lattice. 



Having obtained a simple change of the separated equations (9.15), one can hope that 
there is also a simple modification of the one-dimensional Baxter equations in quantum 
mechanics. Recall that from the works of Sklyanin J3l[ [lj] one knows that the eigenf unctions 
of the quantum Toda lattice Hamiltonian arc given by 

/n— l 
C(\,E)Ml)d\, C{\E) = \{c{\ J1 E). 

3=1 

Here ip\ are renormalized Whittaker functions and the functions c(A, E) satisfy to one- 
dimensional Baxter equation 

P{\) c (A, E) = i n c(X + ih, E) + r™ c (A - ih, E) , 

where P(X) is a trace of the quantum monod rom y matrix T(A). In the classical limit the 



polynomial P(X) enters in the spectral curve (9.6) 



By using a similar approach |3lj, [L4| we can suppose that the eigenfunctions for the dual 
system are expressed in terms of the same Whittaker functions 



i>E(q) = J C(\,E)Ml)d\, 
whereas the corresponding one-dimensional Baxter equation has to be changed 

P(X)c(X,E) =i n (i-E) c(X + ih,E)+r n c(X-ih,E), 
in accordance with the corresponding classical separated equations ( 9.15| ). In the classical 



limit the polynomial P(X) enters in the spectral curve fl9.6[ ). 

Recall that in the Birman-Schwinger formalism we can estimate spectrum and eigen- 
values of the one Hamiltonian H by using known spectrum and eigenvalues of the dual 
Hamiltonian H. So it is interesting to study such a duality in framework of the quantum 
Q-operator theory as an example for the Toda lattice. 

10 The Goryachev-Chaplygin top 

The construction of canonical transformations of the extended phase space proposed for the 
Stackel systems was inspired by the Kepler and Liouvillc results. These transformations 



give rise to the shift of the Lax matrices (6.11) and it allowed us to introduce an integrable 



system dual to the Toda lattice. As above we can try to construct another integrable systems 
starting with known ones by using the mapping of the 2x2 Lax matrices for the Toda lattice 



teD 



We start with any Lax matrix T(X) in the form (3.E). Substituting the known Hamilton 



function H into the mapping (9.11) one calculates a new matrix T and a new Hamilton 
function, which may be tested on integrability. Note that to construct new Lax matrices by 
the rule (6.11 -873) we have to predict the new Hamiltonian H from some external reasons. 



As an example here we consider the Goryachev-Chaplygin top. We introduce coordi- 
nates on the dual space to the Lie algebra e(3) with the standard Lie-Poisson brackets 

\H j ' j j — &ijk ^k 7 \H j 9j t — E-ijk 9k 3 

(10.1) 
{di,9j } = 0, i,j,k= 1,2,3. 
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The orbits on e(3)* are fixed by values of the two Casimir operators C\ — (g,g); C2 — (l,g)- 
The Hamilton function for the Goryachcv-Chaplygin top is equal to 



H=-(lj + ll+4ll)-pl 3 + g 2 , pe 



(10.2) 



It is a completely integrable top on the one-parameter subset of orbits O (C\ — const , Ci = 
0) in e(3)*. The corresponding 2x2 Lax matrix T(A) was obtained by Sklyanin |3C[| . 
According to 1 35 this matrix is closely related to the Lax matrix for the three-particle Toda 
lattice and it may be factored as 



T(A) = Ti(A) T 23 (A) , where T x = 



X 



2Z, 



e q 





and 



T23 — 



A 2 - 2/3A -l\-l\ ie«[A(ffi - i 9 2) - g 3 (h - ih)\ 

ie«[\(gl + ig 2 )-g 3 (h+il2)} 



9l 



(10.3) 



By using the transformation of the Lax matrices (9.11) proposed for the Toda lattices one 
constructs another Lax matrices 



T{\) = T x 



T: 



23 



H -b 





A = v~ 1 A 



which describes a new integrable system on the same one-parameter subset of orbits O. So 



the canonical transformation of the extended phase space (3.4) by v = g 3 , i.e 



H = g\ (H-b), 



(10.4) 



preserves integrability. Moreover, initial and resulting systems are separable in the common 
system of separated variables. By using Maupertuis' principle a similar transformation of 
the extended phase space (10.4) has been obtained in pa]. 



Transformation of the corresponding spectral curves has the expected form 



C: 

C: 



A 2 (g,g) 



= A 3 - p A 2 - 2H A - K , 



X 2 (g,g) + 2H 



= \ 3 -p\ 2 -2b\-K . 



As above, these integrable systems are topologically equivalent. 

Starting from the matrix T 2 3(A) ( |10.3[ ) we can construct 2x2 Lax matrix for the so- 
called Kowalevski-Goryachev-Chaplygin top (see references within pBJ). This Lax matrix is 
closely related to the Lax matrix for the Toda lattice associated with the root system BCi- 
Starting with the induced transformation of the 2x2 Lax matrices for the BC n Toda lattices 
we can construct a new integrable system related to the Kowalevski-Goryachev-Chaplygin 
top. Similar transformations of the extended phase space have been proposed in 1 29 directly 
from Maupertuis' principle. 



11 Conclusion 

The modest aim of this review was to collect some old and new examples of canonical 
transformations of the extended phase space, which map a given integrable system into the 
other integrable system. 

The Sections 4-10 together show that all the examples have many common properties. 
Analysis of these common properties could allow us to join different integrable systems 
into the classes of topologically equivalent systems and to study these classes instead of 



21 



considering of the individual systems. Each of this class of integrable system may be related 
to the class of topologically equivalent n-dimensional Lagrangian submanifolds, which are 
diffeomorphic to the n-dimensional torus. In this approach the different integrable systems 
are associated with the various parametric forms of the common integrable manifold. 

It remains unclear how to construct canonical transformations of the extended phase 
space Me for a given integrable system. Moreover, up to now one does not know all 
consequences of the action of canonical transformations on the Stackel systems, Toda lattices 
and other known integrable systems. 

There is still much to do: to describe modification of bi-hamiltonian structures and 
the Backlund transformations for the Stackel systems and the Toda lattices, to transform 
the corresponding stationary flows of hierarchies of nonlinear evolution equations and to 
consider composed transformations similar to the Kolosoff mapping for the Kowalevski top. 

For all the examples, after canonical transformations of Me, one usually gets dynamical 
r-matrices instead of constant ones. We have to check that these resulting r-matrices satisfy 
to the classical dynamical Yang-Baxter equation and to interpret properly these dynamical 
matrices in a general theory of dynamical r-matrices. 
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